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We identify a new inflationary regime for which more than 60 e-folds are generated classically 
during the waterfall phase occuring after the usual hybrid inflation. By performing a bayesian 
Monte-Carlo-Markov-Chain analysis, this scenario is shown to take place in a large part of the 
parameter space of the model. When this occurs, the observable perturbation modes leave the 
Hubble radius during waterfall inflation. The power spectrum of adiabatic perturbations is red, 
possibly in agreement with CMB constraints. A particular attention has been given to study only 
the regions for which quantum backreactions do not affect the classical dynamics. Implications 
concerning the preheating and the absence of topological defects in our universe are discussed. 
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I. INTRODUCTION 



In the last two decades, precise measurements of Cos- 
mic Microwave Background (CMB) anisotropies have 
provided strong arguments in favor of a phase of ac- 
celerated expansion in the early universe. These argu- 
ments have lead to propose a large number of inflation- 
ary models, most of them based on one or more scalar 
fields slowly rolling along an associated potential, and 
whose predictions on the power spectrum of scalar initial 
perturbations agree with best the measurements of CMB 
anisotropies. Inside the jungle of models, the hybrid class 
is particularly promising because easily embedded in a 
high energy framework like supersymmetry, supergrav- 
ity, or grand unified theories [l[. Contrary to power law 
chaotic models , the energy scale of hybrid inflation can 
be low, and the model do not need super-plankian values 
of the fields. Moreover, it is shown in [3|-Q that initial 
field values in hybrid inflation do not suffer a fine-tuning 
problem, unlike previously thought [H, 0]- Nevertheless, 
the original hybrid model is usually considered as a toy 
model, because its scalar power spectrum calculated in 
the 1-field slow-roll approximation exhibits a slight blue 
tilt, which is disfavored by CMB experiments 0. In this 
paper, we discuss the validity of that calculation. We 
show that in a large part of the parameter space, the 
last 60 e-folds of inflation, relevant for the calculation of 
the observable scalar power spectrum, are actually real- 
ized in a non-trivial way during the waterfall phase, after 
crossing the instability point. 

The original hybrid model of inflation was proposed in 



Refs. [H, Its potential reads 
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The field 4> is the inflaton, ip is the auxiliary Higgs-type 
field, and M, fi, v are three mass parameters. Inflation is 
assumed to be realized in the false- vacuum along the val- 
ley (ip) = 0. In the usual description, inflation ends due 
to a tachyonic instability soon after the inflaton reaches 
a critical value <p c = ^/2v 2 /M. From this point, the clas- 
sical system would evolve quickly toward one of its true 
minima (<p) = 0, and (tp) — ±M whereas in a realistic 
scenario one expects the tachyonic instability to trigger 
a preheating era [TlH3 • 

To determine the scalar power spectrum at the end of 
inflation, it is common usage to restrict the dynamics to 
the effective one-field potential along the valley ip — 



V eB (<l>) = A 4 



(2) 



and to assume that inflation ends abruptly once instabil- 
ity point is reached. Under these hypotheses, the primor- 
dial scalar power spectrum can be easily derived in the 
slow-roll approximiation 19]. It is nearly scale invariant, 
with a spectral tilt 
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where a star means that the quantity is evaluated when 
the pivot mode ko leaves the Hubble radius, that is when 
ko = aH. 

However, let us remind that this formalism is based on 
the belief that inflation ends immediately when the insta- 
bility point <j) c is reached. For original hybrid model and 
its effective potential (Eq. the slow- roll parameters 
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where a prime denotes the derivative with respect to the 
inflaton field. Throughout the paper, we restrict our 
analysis to the case in which inflation occurs in the small 
field phase {<p <C fi). In this regime e\ is extremely small, 
£2 is negative and is the dominant contribution to the 
spectral tilt. Thus the scalar spectral index is generi- 
cally larger than one and the spectrum is blue, which is 
disfavored by WMAP7 observations [20]. 

In this paper, we relax the hypothesis of instantaneous 
end of inflation at the instability point and look for tra- 
jectories performing a phase of inflation during the wa- 
terfall [21]. Tachyonic preheating is not triggered during 
this phase because the effective mass of the adiabatic field 
is small compared to the expansion term. More precisely, 
the exponential growth of perturbation modes, charac- 
teristic of tachyonic preheating, is avoided because the 
Hubble expansion term is dominating the equation gov- 
erning the linear perturbations of ip. In that case, both 
analytical investigations and lattice simulations describ- 
ing the tachyonic preheating do not apply. 

In this paper, the classical dynamics is investigated 
both numerically and analytically by using the adia- 
batic field formalism (22[. The whole potential parame- 
ter space is explored using a Monte-Carlo-Markov-Chains 
(MCMC) method. Regions for which much more than 
60 e-folds are realized after instability are shown to be 
generic. In such cases, observable modes leave the Hub- 
ble radius during waterfall inflation and a modification 
of the predicted scalar spectral index is expected. For 
adiabatic perturbations, the power spectrum is actually 
generically red. 

The potential is very flat near the critical instabil- 
ity point and quantum backreactions could dominate the 



1 Throughout the paper, m p i denotes the physical Planck mass, 
and M p i stands for the reduced Planck mass M p i ~ 0.2m p i ~ 
2.4 x 10 18 GeV. 



classical dynamics. Therefore, a particular attention has 
been given to consider only trajectories not affected by 
quantum stochastic effects. These comprise both the 
quantum backreactions of the adiabatic and the entropic 
transverse field. The classical regime is valid only if clas- 
sical jumps of the adiabatic field are larger than its quan- 
tum fluctuations and if quantum diffusion of the trans- 
verse field [23[ do not increase too much the spread of its 
probability distribution. 

For a short phase of waterfall inflation, some problems 
are put in evidence in [231] . Inflating topological defects 
can induce large-scale perturbations and primordial black 
holes can be formed after inflation. When the water- 
fall phase is much longer, we argue that these problems 
are naturally avoided. Indeed, topological defects are 
so strongly diluted during waterfall inflation that they 
do not affect the observable universe. Primordial black 
holes are expected to form when fractional density per- 
turbations occurring at the phase transition reenter the 
horizon. Thus they affect the observable universe only 
if waterfall inflation lasts less than typically 60 e-folds. 
This is not the case here. 

The paper is organized as follows: section 2 is dedi- 
cated to the dynamics inside the valley, before the in- 
stability point is reached. It is shown that classical os- 
cillations of the waterfall field are quickly dominated by 
its quantum fluctuations. In section 3, we show that 
much more than 60 e-folds can be realized along a water- 
fall trajectory, i.e. after crossing the critical instability 
point. In section 4 the generic character of this effect is 
studied. The dependences on the potential parameters 
and initial conditions are determined by using a MCMC 
method. In the conclusion, important implications for 
hybrid models (e.g. on the formation of topological de- 
fects) are discussed. 



II. FIELD DYNAMICS BEFORE INSTABILITY 

Given an arbitrary set of initial conditions, two clas- 
sical behaviors are possible. Either the trajectory falls 
through one of the global minima of the potential with- 
out inflating. Either it reaches the nearly flat valley along 
(/> = and slow-roll inflation can occur. If the valley 
is reached, trajectories are characterized by damped os- 
cillations in the transverse direction (orthogonal to the 
valley). After some oscillations, slow-roll regime begins 
along the bottom of the valley and a large number of 
e-folds is realised. 

At the critical point of instability 4> c , only a small 
transverse displacement allows inflation to end with a 
waterfall phase. The two competing processes able to 
cause this displacement from the -0 = valley line are the 
remaining classical transverse oscillations and the quan- 
tum fluctuations of the auxiliary field. In this section, it 
is shown that oscillations own generically an amplitude so 
small that they are dominated by quantum fluctuations 
of the auxiliary field. 
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Quantum fluctuations are typically of the order 
Aip qu ~ H/2~k. Primordial nucleosynthesis fixes a rea- 
sonable lower bound on the energy scale of inflation, and 
thus on H through the Friedmann-Lemaitre equations 
(Eq. [6]). On the other hand, measurements of the pri- 
mordial scalar power spectrum amplitude, 



V{k = 0.002/Mpc) ~ 2.43 x 10' 



(5) 



with ei < 0.1, allow to fix a higher bound on H. One 
can thus determine the range of transverse quantum os- 
cillations, 10~ 30 m pl < Atp qu i$ 10" 6 m p i. 

To study the classical oscillations, the two-field dy- 
namics needs to be integrated. For a flat Friedmann- 
Lemaitre-Robertson- Walker metric, the equations gov- 
erning the two-field dynamics are the Friedmann- 
Lemaitre equations, 



decreasing amplitude 
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where C\ and C2 are two integrating constant fixed by 
initial conditions. As an example, it takes about N ~ 40 
c-folds for initial oscillations of amplitude A ~ 10 _3 m p i 
to be reduced by a factor of 10 27 at the minimal level of 
quantum fluctuations. 

Given this time scale, the assumption that </> is con- 
stant can be justified a posteriori. Indeed, in the slow-roll 
approximation, straightforward manipulations [§| give 
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as well as the Klein- Gordon equations 

<p + 3H(t) H — — = , 

ocp 
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ip + oHip H — = . 



(7) 



where H = a /a is the Hubble parameter, a is the scale 
factor and a dot denotes derivative with respect to cosmic 
time. 

In the regime of small classical oscillations -C M, 
at inflaton values sufficiently larger than the critical one 
but still in the small field phase, that is <fi c <C <fi <C /z, the 
potential is well approximated by 



V((j),ip) ~ A 4 
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We can assume that the inflaton field is slow-rolling 
along the valley "0 = 0, such that 



H 



1 



V3M, 



pi 



VA 4 (l + 2 /^ 2 ) 



(9) 



In a short time scale, the inflaton </>, and thus the Hub- 
ble parameter H, can be assumed to be constant. As a 
consequence, the Klein-Gordon equation for the auxiliary 
field now read 

i> + 3f— L-,/A*(l + <W) + = . (10) 

It has a simple oscillating solution with exponentially 
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where (f>\ is the initial inflaton value. Therefore, if <C 
/j,, classical oscillations become dominated by quantum 
fluctuations in a range of inflaton value 
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which is typically very small for a nearly flat valley. Thus 
the classical oscillations of ip are expected to be domi- 
nated by quantum fluctuations after a very small range 
of variation for <p. 

To go beyond this approximation and determine, for 
the full potential, how generic are trajectories whose clas- 
sical oscillations become dominated by transverse quan- 
tum fluctuations, the classical 2-field dynamics have been 
integrated numerically. We have followed the method 
used in Q and run an identical Monte-Carlo-Markov- 
Chains analysis of the 7D space of initial field values, 
initial velocities and potential parameters. The result of 
this analysis is that around 99.9% of trajectories trapped 
inside the valley perform transverse oscillations whose 
amplitude is below the most restrictive level of transverse 
quantum fluctuations A-0 qu ~ 10 _30 to p i. 

From these considerations, at instability, the slight 
transverse displacement essential for a waterfall phase 
to take place is not supplied by classical oscillations of 
the auxiliary field but by its quantum fluctuations. In 
the following sections, the waterfall phase will be studied 
classically taking initial values </>; = 4> c and ipi — Atp qu . 
We will assume initial field velocities given by the slow- 
roll approximation. Actually, due to the slow-roll attrac- 
tor, different choices of initial velocities only marginally 
influence the resulting waterfall dynamics. 
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III. HOW HYBRID INFLATION ENDS 

Before to study the waterfall phase, it must be veri- 
fied that the classical dynamics is valid and not spoiled 
by quantum backreactions of both the adiabatic and the 
entropic fields. 



A. Quantum backreactions 

The collective evolution of the fields can be described 
by the adiabatic field a, defined hereafter in Eq. (|36|) . It 
is very light and its classical evolution is valid if classical 
jumps are larger than the quantum fluctuation scale, that 
is 



<* A H 

A(T r i = — > Aer ml ~ — 



which is equivalent to 



ei(cr) > r 

7rm pl 



(14) 



(15) 



Thus we pay a particular attention to only consider wa- 
terfall trajectories along which this condition remains 
true. 

At the critical instability point, the classical value of 
the transverse field is about 0, and the transverse quan- 
tum fluctuations will determine on which side the system 
will evolve towards. The overall dynamics remains clas- 
sical due to the (f> field evolution. However, it must be 
ensured that the quantum backreactions of the transverse 
field do not push the field evolution far from the valley 
line ip = 0. Such effects would modify strongly the dy- 
namics such that the waterfall phase would take place 
in a low number of e- folds. In other words, it must be 
ensured that the spread of the probability distribution of 
the auxiliary field does not becomes much larger than its 
classical value during the waterfall. 

The coarsed-grained auxiliary field is described by a 
Klein-Gordon equation in which a random noise field £(i) 
is added. This term acts as a classical stochastic source 
term. In the slow-roll approximation, the evolution is 
given by the first order Langevin equation 
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The two-points correlation function of the noise field 
obeys to 

<£(«)> = 0, mm) = 5(t-t>) . as) 



When the dynamics is dominated by the <j> evolution, 
this equation can be integrated exactly. Under a change 
of variable [23| . x = exp [— 2r(N — N c )], where 



V 



(19) 



and where N c is the number of e-folds at critical point of 
instability, one has 



dtp 
dx 
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The exact solution to this equation is 

ip(x) = C exp [C2X — C% In x] 
— C\ exp [C2X — C2 In x] 



(20) 



J exp [-C 2 x' + C 2 In a/] ^(x')dx' , 



(21) 



where C\ = i? 1/2 /47rr, C2 = 2/M 2 r and C is a constant 
of integration. Taking the two point correlation function 
and assuming an initial delta distribution for -0 at <f> ^> 
4> Cl one can obtain 



H 2 



exp(x) 



ax 



r(a, ax) , 



(22) 



where a = AM 2 l /M 2 r and T is the upper incomplete 
gamma function. 

In the following, we will consider large values of fi 
and relatively small values of M compared to the Planck 
mass, such that r ~ 2M 2 J [± 2 and a ~ 2fi 2 /M 2 > 1. At 
instability, x — 1 and one thus has 
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By using recurrence relations as well as the asymptotic 
behavior of the T function, one can find 



/ e \ u 

(-J xT(u,u) 



IT 1 
2J^L 



when u — > 00 , (24) 



such that 



(^(x = l)) 
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32n 3 / 2 M 2 l 



(25) 



For instance, for the parameter values of FigfTJ fj, = 
636.4 m p i,M = 0.03 m p i, one obtains v/ (ip 2 ) ~ 2H It 
will be shown later in the paper (see Sec lIV| that when 
the tachyonic preheating is triggered and forces inflation 
to end, one has x < 1 such that (ip 2 (x en <i)} ~ {4> 2 {x = 1)) 
~ H . Therefore, after the critical instability, all the time 
the field dynamics is governed mainly by the <fi evolution, 
the standard deviation of the transverse field distribution 
around its classical value remains at the same order. 
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Notice that an identical result can be obtained using 
the formalism developped in [2^ - [28| . Since the auxiliary 
field is well anchored at its minimum tp = before the 
waterfall, it can be regarded as the same as its fluctua- 
tion, tp — Sip. As for the stochastic formalism, one can 
assume that the <p field evolves independently according 
to Eq. (|12[) . By Fourier expanding Sip and neglecting non 
linear terms, one obtains the mode evolution equation 



a 2 H 2 



12 



M, 



pi 
M 2 



-2r(N-N c ) 



Stp k = o 

(26) 

where a prime denotes derivative with respect to the 
number of e- folds. Following Ref. [24j . in the high fre- 
quency limit this equation can be solved in terms of the 
WKB approximation. In the low frequency limit its ex- 
act solution is a combination of the Hankel functions 
of first and second kind. Let introduce k c = a c H the 
mode leaving the Hubble radius at the critical instabil- 
ity, and n = N — N c . Then, under the assumption that 
12Afpj/Af 2 ^ 1, one finds for small scales modes near 
the instability k 3> fc c exp(n), 



\6Mk,n)\ = 



H 



x exp [ 



3/2 
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and for large scales modes k <§; k c exp(n), 
\SMk,n)\= 11 



v /2afc3 
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where A = 3 2 / 3 r(2/3)a- 1 / 6 /(20F) is a typically order 
unity factor, and 
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In the regime 2rn <C 1, the modes which become tachy- 
onic satisfy 



< a 2 n e 2n 



(30) 



In (24|, it is then assumed that a 1 and n ~ O(l) 
to find that the quantum back-reactions from the small 
scales entropy perturbations dominate and force infla- 
tion to end quickly after waterfall instability. We will be 
here interested by the opposite case, a < 1. As shown 
later in Sec lIVi the total number of e-folds that can be 
realized classically between the instability and the begin- 
ning of the tachyonic preheating is larger than 60 and is 
roughly given by n ~ ^l 2 M 2 . Therefore the tachyonic 
modes stand super-Hubble during all this phase. In that 



case, the variance of Sip is dominated by the large scale 
mode contribution 



d 3 fc 
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One obtains just after the critical instability 
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which is identical to Eq. (1251) up to an order unity nu- 
merical factor 2 . 

In the section UlI CI initial values of tp at critical insta- 
bility are taken to follow a gaussian random distribution 
that verify Eq. (|23|) . From this point, the classical value 
of tp moves away its initial amplitude and increases such 
that it becomes quickly much larger than its quantum 
fluctuations, even if the overall dynamics is still governed 
mainly by the <p evolution. Therefore, the classical dy- 
namics is not spoiled by transverse quantum fluctuations. 



B. Transverse field gradient contribution 

Another effect susceptible to spoil the classical homo- 
geneous dynamics is the backreaction due to the trans- 
verse field gradient contribution to the energy density. 

Assuming the statistical homogeneity, the mean-square 
value of transverse field gradient after smoothing on a 
length L = 1 /aH is given by 



1 



(27T) 3 
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During the waterfall, tpk is given by Eq. ([28 
gration over the modes, one obtains 



(33) 



After inte- 



H 



(34) 



Since (\tp\ 2 ) ~ H 2 during the waterfall in the regime of 
interest, the transverse field gradient contribution to the 
energy density is negligible compared to the potential 
term V ~ 3H 2 /M 2 V The background dynamics thus 
remains mostly homogeneous. 



C. The classical regime 

Once the instability point is reached, trajectories de- 
viates from the valley line and fall through one of the 



2 Notice that a similar result can be obtained for n = 0, from the 
Eq. (2.27) and (2.30) of Ref. In that case, an additional 

factor a -1 / 6 is obtained, but it can be due to the matching 
problem between large and small scale modes. 
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global minima (<f> = 0, ip = ±M). Following [15|, we can 
assume that the auxiliary field reacts faster than the in- 
flaton field such that trajectories follow the ellipse defined 
by the minima in the ip direction 



dip 



= with 



< 6 < 



M 2 



= 1 



(35) 

The small-field type effective potential defined by this 
ellipse is nearly flat around the instability point, where its 
curvature is negative. Thus a phase of inflation should be 
possible near the top. The adiabatic field a, introduced 
in [22|, defined such that 



a = y (p 2 + tp 2 , 



(36) 



describes the collective evolution of the fields along the 
classical trajectory. One can determine its equation of 
motion 



a + 3H& + V a = 



where 



<pdV ipdV 



a d(p a dip 
On the ellipse of Eq. (|35|) . one obtains 



(37) 



(38) 



V„ = A 
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(39) 



where <f> is related to the adiabatic field through the re- 
lation 



a(4>) 




(40) 



Like for an effective f-field model, the slow-roll regime 
can be assumed and the slow-roll parameters can be in- 
troduced. 

Let remark again that this approach is valid under the 
assumption that trajectories follow the above defined el- 
lipse. In practice, at the critical instability point, the 
gradient of the potential is along the <j> direction. The 
field evolution first follow this direction and thus does 
not follow exactly the ellipse of minima. Therefore the 
predictions are expected to be modified more or less im- 
portantly. These modifications are studied by solving 
numerically the exact classical dynamics. 

Fig. [T] shows a typical trajectory, for typical poten- 
tial parameters. Starting integration at instability with 
ipi = 10 _12 m p i ensures that the dynamics of the adia- 
batic field is not dominated by its quantum stochastic 
effects. Indeed, the initial transverse displacement cor- 
responds to the scale of quantum transverse fluctuations 
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FIG. 1: Top: typical trajectory (solid line) in the {cj>,ip) 
plane {(j> and tp are in m p i units), for initial field values 
(pi = 4>c = 0-03m p i,^j = 10 -12 m p i and potential parameters 
M — 0.03m p i,/i = 636.4m p i. The ellipse of minima of Eg . (|35|) 
is also represented (dashed line). From left to right, the points 
indicate where ei = f _3 /10 _2 /0.f /I . Center: zoom around 
4> c for the same trajectory. From left to right, the points in- 
dicate where n B = 1./0. 97/0.91/0. 65. Bottom: Number of 
e-folds realised along this trajectory, from the critical insta- 
bility point. 



given by Eq. (|23|) . Thus for our set of parameters, one 
can read 



<t>lm 2 vl 
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(41) 
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and thus quantum jumps of the adiabatic field are much 
smaller than classical evolution during one Hubble time. 

More than 600 e-folds are found to be realised before 
inflation ends. Therefore, 60 e-folds before the effective 
end of inflation, when observable modes cross the Hub- 
ble radius, instability has already been crossed and the 
trajectory follow a non-trivial evolution. 

The effective potential of Eq. (|39l) is thus an ideal case. 
Two regimes during which inflation is possible are iden- 
tified numerically: 

1. PHASE I: At <j> = <p c , it is assumed that each tra- 
jectory is slightly displaced in the i\> direction. Tra- 
jectories first exhibit a non trivial behavior. Driven 
by their velocity along 0, they first follow the slope 
in the <p direction before turning, following roughly 
the gradient of the potential until the ellipse de- 
fined in Eq. (|35|) is reached. As shown in FigfT] a 
large number of e-folds can be realized during this 
first phase. 

2. PHASE II: Trajectories reach and mostly follow the 
above defined ellipse. A large number of e-folds is 
realized if the effective potential along the ellipse is 
sufficiently flat. 

Finally, if in the standard approach the tachyonic pre- 
heating begins at phase transition, in the case of waterfall 
inflation it can only occurs at the end of inflation, simi- 
larly to what happens for the new inflation model |29| . 

The full calculation of the primordial power spectrum, 
including entropic perturbations, is beyong the scope of 
this paper. However, along the ellipse, the scalar spectral 
index for adiabatic modes can be evaluated in the slow- 
roll approximation. One has to evaluate the field value 
at Hubble exit of observable modes, that is when 



N{(j>) 



Then the spectral index is directly determined with 
Eq. (U). It is generically red. The spectral index of the 
adiabatic power spectrum can also be determined for nu- 
merical trajectories. It has been plotted along the tra- 
jectory of Fig. [T] as well as for a grid in the parameter 
space (fx,M), see Fig.[S] 

In this section, some waterfall trajectories leading to 
more than 60-folds have been found to exist. But before 
to draw conclusions, it is essential to measure how generic 
such trajectories are in the parameter space. This is the 
point of the following section, in which the full poten- 
tial parameter space will be explored using a statistical 
MCMC method. 




through its three parameters M, fx, v. In order to draw 
more easily the physical interpretations, we have replaced 
the parameter v by the position of the instability point 
<j) c . The dynamics depends also on the initial value of 
the auxiliary field. It is given by the scale of its quantum 
fluctuations given by Eq. (|23[) at instability point when 
it acquires a large mass. So it is related to the potential 
parameter A through the Friedmann-Lemaitre equation. 

To explore this 4D space, we have used a Monte-Carlo- 
Markov- Chains method. Flat priors have been chosen on 
the logarithm of these parameters, in order to not favour 
any precise scale. The chosen ranges of parameters are 
the following: 



0.3m p i < fx < I0 4 m p i 
I0" 6 to p i < M < Mj 



'pi 

I0~ 6 to p i < <j) c < M p i 
10- 60 m 4 



pi 



< A 4 < 10 



-12 4 



(43) 
(44) 
(45) 
(46) 



The lower bound on /i comes from its posterior probabil- 
ity distribution [3] to generate sufficiently long inflation 
inside the valley from arbitrary subplankian initial con- 
ditions. Upper bounds on M and </> c stand because we 
only consider the dynamics at field values smaller than 
the reduced Planck mass 3 . Lower bounds on M and 4>c 
and upper bound on /i are arbitrary for numerical conve- 
nience. Prior bounds on the parameter A are such that 
the present constraints on the energy scale of inflation, 
given by nucleosynthesis and observations of the CMB, 
are respected. 

Trajectories are integrated from the instability point. 
Initial inflaton value is thus <j){ — <p c . Initial auxiliary 
field values are assumed to follow a gaussian distribution 
around ip = 0, with a dispersion given by Eq. (|23|) . In 
order to avoid strong quantum backreactions of the adi- 
abatic field, a hard prior coming from Eq. [IS] is enforced, 



> 



H 2 



A 4 



(47) 



'pi 



such that each trajectory that do not verify this condi- 
tion is excluded of the Markov chain. Trajectories are 
integrated until e\ — 1, the end of inflation. The accep- 
tance condition to a Markov chain is a realization of at 
least 60 e-folds after </> c . 

Marginalized posterior probability density distribu- 
tions, normalized such that the area under each distri- 
bution is 1, are shown in Figs. [2]and|4j Posterior dis- 
tributions of the parameters M and \i are degenerated, 
but the posterior distribution of their product fj,M 4 (see 



IV. EXPLORATION OF PARAMETER SPACE 

The number of e-folds generated after crossing the 
instability point depends on the form of the potential 



3 this restriction appears to be natural in SUSY inspired models, 
for which the potential is lift up at super-Planckian field values 
due to radiative corrections. 

4 in that case, the parameter fj, is replaced by the product Mfi in 
the MCMC simulation, using a flat prior on the log of fiM. 
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FIG. 2: Marginalized posterior probability distributions of 
the potential parameters fj, (top) and M (bottom). The ver- 
tical axis is normalized such that the total area under the 
distribution is 1. 



FIG. 4: Marginalized posterior probability distributions of 
the critical point of instability <j} c (top) and the normalizing 
parameter A 4 (bottom). 




FIG. 3: Marginalized posterior probability distribution of the 
product Mfi. 



FigEJ) is prior independent. A bound on this combination 
is: 



log 



>0.21 95%C.L. 



(48) 



This bound can be rewritten M/m p \ > m p i//z and ex- 
plained intuitively. A large number of e-folds have to be 
generated after the instability point, before the magni- 



tude of the effective negative mass of the auxiliary field 

-A 2 / 



M 



(49) 



increases and becomes larger than H, forcing inflation to 
stop. Following (23[, this happens in the range 



> <p > 



'1 - 



M 2 



(50) 



'pi 



During this period, the slow-roll approximation for <p is 
valid and the number of e-folds generated is given by 
Eq. ([l"2j) . In the limit <j) <C using the range of <f> ob- 
tained in Eq. (|50l) . straightforward manipulations give 
AN ~ n 2 M 2 /4M^ and thus the number of e-folds is 
roughly fixed by the combination of the parameters fj, 
and M. From this reasoning comes also the argument 
that, at the end of the Phase I, x ~ exp[-M 2 /M^j] - 1 . 

Since the analyis is restricted to the sub-planckian field 
dynamics, and thus to sub-planckian values of M, there 
exist a gap between 



0.3m p i < /i < 6r 



ipl, 



(51) 



for which the number of e-folds generated is less than 60. 
In this regime, the slope of the potential in the <f> direction 



9 



is too large to generate a sufficient number of e-folds in 
Phase I, before the negative mass of the auxiliary field 
pushes the trajectory away from the -0 = line. 

The posterior probability distribution of <p c is nearly 
flat and this parameter does not influence significatively 
the duration of the waterfall inflationary phase. 

The marginalized posterior probability distribution of 
A 4 is almost flat and decreases for values corresponding 
to the highest energy scales of inflation, without becom- 
ing negligible. This suppression is not only due to the 
hard prior. More important is the effect of larger initial 
values of the auxiliary field on the dynamics. When clas- 
sical trajectories are initially far away from the tp — 
line, the phase I is reduced and less efficient to generate 
a large number of e-folds. 

The MCMC analysis therefore provides an explicit an- 
swer to the question of how generic are trajectories realiz- 
ing a large number of e-folds after the critical instability 
point. They are found to occupy a large part of the pa- 
rameter space, gathered in the region where Eq. (|48|) is 
verified. Since A 4 is directly linked to the energy scale 
of inflation, waterfall inflation is also found to be more 
favorable at low energy compared to the Planck scale. 

Finally, let remark that the posterior probabil- 
ity distributions obtained with the MCMC method 
should be combined to the probability distributions of 
4>i, ipi, <f>i,"4>i, M, fi, 4> c , obtained in [3], related to the prob- 
ability for trajectories initially outside the valley to reach 
the slow-roll attractor inside the inflationary valley. 

All these results stand for the original hybrid model, 
but the general features are expected to be reproducible 
with more or less efficiency, for all models in which in- 
flation mostly occur in a nearly flat valley and end due 
to a tachyonic instability, like in many SUSY realizations 
(e.g. F-term hybrid model [30l|). 

Notice that when Eq. 25] is not verified, then the stan- 
dart mechanism do work: namely inflation stops soon 
after <j> c . 



V. CONCLUSION AND DISCUSSION 

In hybrid models, the usual description assumes slow- 
roll inflation along a nearly flat valley ending quasi in- 
stantaneously due to a tachyonic instability triggering 
the tachyonic preheating [llTfl8l |. 

In this paper, the waterfall phase has been studied in a 
regime during which inflation continues for a long time. 
It has been shown that much more than 60 e-folds can 
be realised classically during the waterfall, after cross- 
ing the instability point. Particular attention has been 
given to study regions in the parameter space where the 
classical dynamics is valid and not spoiled by quantum 
backreactions of adiabatic and entropic fields. 

Observable modes leave the Hubble radius when the 
effective potential is very flat with a negative curvature. 
Instead of blue, the power spectrum of adiabatic pertur- 
bations is thus generically red. However, the calculation 




FIG. 5: 200 x 200 grid of spectral index values of the power 
spectrum of adiabatic perturbations, in the plane (/x, M) for 
the exact classical dynamics, with <j) c = 10 -3 m p i. Black re- 
gion correspond to trajectories leading to no more than 60 
e-folds after instability. 



does not include the contribution of potentially observ- 
able iso-curvature modes. This contribution was shown 
to be negligible in recent developments [24H28j |. but in 
these references a fast waterfall phase is assumed and 
their conclusion does not apply to the regime of poten- 
tial parameters studied here. The full numerical calcula- 
tion of the primordial power spectrum should be realised 
soon in [31J. The effect on the non-gaussianities pro- 
duced during the tachyonic preheating [32|, l33| could be 
also important. 

Therefore, it is premature to conclude that original hy- 
brid model is strongly disfavored by CMB experiments. 
All the more since a bayesian MCMC analysis demon- 
strates that such trajectories are generic in a large part 
of the potential parameters space. 

Our result is expected to be present in all hybrid mod- 
els with an inflationary valley and a tachyonic instability. 
In particular, a preliminary study of F-term supersym- 
metric hybrid model ^30] indicates that it exhibits a sim- 
ilar dynamics in some part of its parameter space. 

These observations may have also an important im- 
pact on questions related to the end of inflation. In par- 
ticular, if a large number of e-folds occur after symme- 
try breaking, the eventually formed topological defects 
will be diluted by expansion and thus will not affect our 
observable universe. Therefore, some works constrain- 
ing the schemes of symmetry breaking in grand unified 
theories with topological defects [34[ may be reviewed. 
Our results should have also some impact on tachyonic 
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preheating. In our knowledge, all the previous studies 
of tachyonic preheating in hybrid inflation [13l - ll5j ne- 
glect expansion. Lattice simulations indicates that cos- 
mic strings and domain walls strongly affect the way pre- 
heating phase occur. However if such defects are diluted 
by a phase of inflation after symmetry breaking, lattice 
simulations should be updated to include expansion or to 
start after waterfall inflation, like in new inflation models 

Finally, let comment about stochastic effects. For 1- 
field effective potential of hybrid inflation, these were 
found to not affect the classical trajectories along the 
valley [H, Hfil- Authors also notice that in small field 
inflation, a stable solution of eternal inflation should ex- 
ist at the top of the potential. At cf> ~ 0, the hybrid 
potential along ip reduces to a small field type and thus 



this observation should apply to the waterfall phase in 
hybrid inflation. In particular, it would be interesting to 
determine, in a full 2-field approach of stochastic effects, 
how the field dynamics is affected if initial conditions for 
the waterfall are taken along the border of the stochastic 
patch. 
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